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ABSTRACT 

Relativistic coUisionless shocks are believed to be efficient particle accelerators. Nonlinear outcome 
of the interaction of accelerated particles that run ahead of the shock, the so-called "precursor" , with 
the unperturbed plasma of the shock upstream, is thought to facilitate additional acceleration of 
these particles and to possibly modify the hydrodynamic structure of the shock. We explore here 
the linear growth of kinetic modes appearing in the precursor-upstream interaction in relativistic 
shocks propagating in non and weakly magnetized plasmas: electrostatic two-stream parallel mode 
and electrostatic oblique modes. The physics of the parallel and oblique modes is similar, and thus, 
we refer to the entire spectrum of electrostatic modes as "two-stream-like." These modes are of 
particular interest because they are the fastest growing modes known in this type of system. Using 
a simplified distribution function for a dilute ultra-relativistic beam that is relativistically hot in 
its own rest frame, yet has momenta that are narrowly coUimated in the frame of the cold upstream 
plasma into which it propagates, we identify the fastest growing mode in the full /c-space and calculate 
its growth rate. We consider all types of plasma (pairs and ions-electrons) and beam (charged and 
charge-neutral). We find that unstable electrostatic modes are present in any type of plasma and 
for any shock parameters. We further find that two modes, one parallel {k± = 0) and the other one 
oblique {k± ~ A;||), are competing for dominance and that either one may dominate the growth rate in 
different regions of the phase space. The dominant mode is determined mostly by the perpendicular 
spread of the accelerated particle momenta in the upstream frame, which reflects the shock Lorentz 
factor. The parallel mode becomes more dominant in shocks with lower Lorentz factors (i.e., with 
larger momentum spreads). We briefly discuss possible implications of our results for external shocks 
in gamma-ray burst sources. 

Subject headings: acceleration of particles — cosmic rays — plasma — shock waves 
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1. INTRODUCTION 

Cosmic rays are believed to be acce lerated in astrophysic al collisionless shocks. The leading accelerati on process is dif- 
fusive shock acceleration (DSA; e.g., Beir'1978; Blandford fc OstrikeilUQTStlBlandford k EichledflQSTi ) . where charged 
particles are accelerated by crossing the shock back and forth. Thus, an integral element of this picture is a precursor of 
accelerated particles that runs ahead of the shock, and interacts with the incoming, and mostly unperturbed, upstream 
plasma. The upstream plasma must interact with the precursor in such a way as to deflect the accelerated particles 
back into the shock and thus close the Fermi cycle and enable farther acceleration. It also provides an early pertur- 
bation to the upstream plasma that may, through a transfer of energy and momentum, and throu gh an exc itation of 
electromagnetic perturbations in the plasma, modify the shock structure (e.g.. jBlandford 1980: D r urv fc Voelk 19811. 
If a different acceleration mechanism is at work (e.g., shock-drift acceleration: I Webb! Axford fc Terasawal 19830 . then 
the accelerated particles should similarly cross the shock into the upstream, thus producing a precursor, but these 
particles are not necessarily deflected back into the shock. 

The general picture described above should apply to both relativistic and non-relativistic shocks and in various types 
of upstream plasma and precursors, independelty of the plasma magnetization level and composition. The details of 
the precursor interaction, however, will depend on the specific system properties. Here we explore a family of such 
systems that are thought to occur in astrophysical environments: relativistic collisionless shocks that propagate into 
unmagnetized or weakly magnetized plasma of various compositions. Until recently, theoretical understanding of 
particle acceleration in such shoc ks was based mostly on analytic or semianalytic prescriptions (e.g.. lAchterberg et al. I 
120011 : [Ellison fc Dou ble 2002; Ni emiec fc Ostrowskill2006l and references therein), in which the structure of the medium 
that scatters the accelerated particles was assumed ad hoc. A major breakthrough was achieved recently, when two- 
dimensional first-principles particle-in-cell (PIC) simulations of unmagnetized pair and ion-elect ron shock s provided 
direct evidence for the acceleration of particles via DSA and for the formation of a shock precursor (|Spitkovsk v 2008a^lJ; 
iKeshet et al. II2009I : iMartins et al. |[2009ll . As expected, the numerical simulations show that a precursor of accelerated 
particles runs ahead of the shock and interacts, by exciting collective modes, with the initially undisturbed upstream 
plasma. While opening a new window into the physical processes that take place in these shocks, PIC simulations 
are still rather limited in size and duration and an analytic treatment of the processes that shape the shock is still 
required to extrapolate between the simulations and real astrophysical systems. Moreover, existing simulations have 
not converged to a steady state shock structure. The maximum energy of the accelerated particles continues to grow 
with the duration of the simulation and so does the spatial extent of the precursor, as measured from the shock 
transition. The properties of the leading edge of the precursor also keep evolving in the simulations, e.g., the particle 
density at the precursor's leading edge keeps dropping as the latter progresses away from the shock transition. This 
implies that the kinetic instabilities that will develop at the leading edge, on time and length scales that are not 
sampled by existing simulations, may be very different than those that are observed n ow in exist ing simulations. 

There are several known instabilities that operate in interpenetrating plasmas (see iBretl 120091 for an overview and 
references). Three specific types of instabilities excited in shock precursors could be important in relativistic unmag- 
netized or weakly magnetized shocks. These are the kinetic two-stream-like instability, the ki netic f ilamen tation (or 
transverse Weibel) instability, and the magnetohydrodynamic streaming instability identified bv I Belli ()2004[ ). The Bell 
instability develops only if the upstream plasma has some pre-existing seed magnetic field. The filamentation instabil- 
ity is suppressed if the accelerated particle precursor is dilute (i.e., its particle density is sufficiently small compared 
to the upstream p lasma density), and its perpendicular momentum spread in the upstream frame is large enough 
(jSilva et al. 1 120021 ). The two-stream-like instability, which we define to include the classical two-stream instability 
as well as its oblique variant, is the only one that grows under any circumstances, independent of the density and 
momentum spread in the upstream frame of the shock precursor. Because of the two-stream-like instability's universal 
occurrence, we dedicate this work to identifying its fastest growing modes and discussing its potential role in relativistic 
unmagnetized shocks in astrophysical and numerical settings. 

The two-stream-like and filan ientation instabilities (sometimes referred to as "Weibel" instabilitv lWeibeilll959l : iFriedl 
[Tm LMoiseev & SagdeevHWeS") arise f rom t he same branch of the dispersion relation and compete with each other 
(e.g.,|Brct 2009; Michno fc Schlickciseil l2010[ ). The two-stream parallel unstable mode is found when the wave vector 
is aligned with the flow, while the filamentation mode has its wave vector perpendicular to the flow. For intermediate 
orientations, the so-called "oblique" modes are equally unstable. Since no magnetic field is generated, these oblique 
modes are electrostatic, and are very similar to the parallel two-stream mode, explaining why they can be labeled 
"two-stream-like." For Maxwellian electron beam/plasma system in the non-relativistic regime, two-stream modes 
grow faster than any other unstable mode. In rather high density relativistic beams, filamentation modes tend to 
dominate the linear phase. Oblique modes be come importa nt in the dilute r elativistic beam regime that is the focus 
of the present work (jFamberg. Shapiro fc Shev chcnko 1970: lBret et al. II2008D . 

Motivated by astrophysical considerations, we consider a relativistic beam that is seen in the rest frame of the 
background plasma to be having particle Lorentz factors narrowly distributed around 7, and momenta that are directed 
within an angle of ~ F"^ from the beam axis, where F is the shock Lorentz factor. This is an ultra- relativistic beam 
that is re lativistically hot in its own rest frame, namely 7 3> F ^ 1. Unstable modes in this regime were explored very 
recently (jLemoine fc Pellet"ieill2011|) . but only for modes that satisfy fc||C = Wpi//? , where Wpi is the background plasma 
frequency and /3 is the velocity of beam particles. Here we explore the entire (fc|| , k±) space, including the pure parallel 
modes k± — 0, and derive the fastest growing mode and its growth rate under the assumption that the background 
plasma remains cold. Our results are applicable to all types of unmagnetized shocks and precursors, including beams 
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containing ions and electrons, ions only, electrons only, or pairs, and upstream plasmas containing ions and electrons 
or pairs. 

We do not go in the present work beyond the linear regime. The nonlinear phase should be explored in the future 
(possibly with numerical simulations) , since it is on the detailed character of the nonlinear development that a possible 
significant heating by, or momentum transfer from, the accelerated particles, is contingent. Similarly, the nonlinear 
development determines whether the resulting electrostatic field can contribute accelerated particle scattering in the 
shock upstream. 

The paper is organized as follows. In §2 we discuss the astrophysical context. In §3.1 we define our model distribution 
function and derive the growth rates of the fastest growing modes. In §3.2 we provide a simple physical description 
of the instability. In §4 we discuss possible implications of our results for GRB external shocks. In §5 we provide our 
conclusions. 

2. THE ASTROPHYSICAL SETTING 

Relativistic coUisionless shocks are thought to occur in various astrophysical environments, including in GRBs, pulsar 
wind nebulae, active galactic nuclei, and microquasars. While the type of plasma (ions and electrons or pairs) and 
the magnetization level is debated in most of these systems, it is likely that astrophysical relativistic shocks take place 
in both plasma types and over a large range of magnetization levels. Among these environments, one of the systems 
where the plasma properties are strongly constrained by observations, is the external shock in GRBs. Here an ultra- 
relativistic blast wave starts with very high Lorentz factors (up to F ~ 1000 and maybe even more) and decelerates 
with time to N ewtoni an velocities. It propagates into a very weakly magnetized circum-burst proton-electron plasma 
(for reviews see iPiran I [2004: Mcszaros 2006; Nakar 2007). Another environment were the plasma type is constrained is 
the termination shock of pulsar wind nebulae, which is driven into a pair plasma (although s ome ions may be present), 
but i n this case, the magnetization level of the plasma at the shock is debated (for review see lKirk. Lvubarskv fc Petril 
EOQi). 

Here, we consider coUisionless relativistic shocks of both plasma types, which are non-magnetized or only weakly 
magnetized. The accelerated particles crossing into the shock upstream constitute a beam of relativistic ions and/or 
electrons, or relativistic pairs. The acceleration produces a spectrum of energies that is typically approximated with 
a power-law, where the largest number of the accelerated particles, and most of the beam energy contained in the 
accelerated particles, resides at the low-energy end of the spectrum. Assuming that the particles are accelerated by 
DSA, the ones that run ahead of the shock must be deflected by interacting with the upstream plasma. Lower energy 
particles are deflected more easily than those with higher energies. Therefore the particles with higher energies can 
travel farther ahead of the shock and at any distance from the shock the accelerated beam is dominated, both in 
particle number and total energy, by the lowest energy particles that reach that distance. These are also the particles 
that we expect to be most susceptible to kinetic instabilities and are thus, those that will dominate the growth rate 
of the two-stream-like instability^ Therefore, we approximate the beam at a given distance from the shock as being 
monoenergetic with particles Lorentz factor 7, as observed in the upstream frame, which is the lowest Lorentz factor 
of beam particles at that distance. The rest frame of the beam, i.e., the frame in which the beam is isotropic, moves 
with a Lorentz factor ~ T with respect to the upstream. The beam particles are also highly relativistic in the beam 
rest frame, implying 7/r ^ 1. Thus, we work in the regime in which 7 3> L ^ 1. 

The charge neutrality of the beam depends on the specific setting. A beam of pairs will always be charge-neutral, but 
a proton-electron beam may be either charge-neutral or positively c harged. If electrons and protons are accelerated 
with similar efficiencies (as numerical simulations suggest, see, e.g., lSpitkovs"kvl l2008al) and radiative losses can be 
neglected, then at any distance in front of the shock the beam is charge-neutral and contains particles with similar 
energies (i.e., rrie^fe rnpjp). Therefore, since in the limit 7e,7p S> 1 the beam behavior depends only on the product 
7717, the beam is similar to a beam of pairs. If, on the other hand, electrons are not accelerated as efficiently as protons, 
either due to less efficient injection into the DSA cycle, or due to higher radiative losses, the accelerated particle beam 
is positively charged. Note that the electrons and protons of an initially cold upstream always behave differently. 

An instability can have a strong impact on the structure of an astrophysical system if its amplitude becomes 
nonlinear. In order to do so it must have enough time to grow, namely, the time that the instability has to grow, igrow, 
must be significantly longer than the instability growth time, given by the inverse growth rate 1/S. In the context of 
our setup, the value of tgrow for an upstream fluid element is the time between the first encounter of the fluid element 
with the accelerated particle precursor and the time that this element is swept up by the shock transition. If the shock 
is spherical and its radius is R, then tgrow ^ R/cT'^, where tgrow is measured in the upstream rest frame. In a planar 
geometry, which is typically used in numerical simulations, igrow < V^^, where t is the age of the shock (both times 
are measured in the upstream frame). Note that a beam particle spends a time much longer than tgiow in the shock 
upstream. In spherical geometry it can be as long as ~ R/c and in planar geometry it can be as long as t. Therefore, 
beam particles have much more time to be affected by the interaction with the upstream. 

3. THE FASTEST GROWING MODE 
3.1. The dispersion relation 

^ Our analysis validates this expectation. Below, we find that the instability growth rate is slower for more energetic particles (higher 7) 
and for beams with lower densities. This result supports our monoenergetic approximation as it implies that particles with much higher 
Lorentz factor than 7 are reacting to the the background plasma on much longer time scales in the linear regime and can therefore be 
neglected. 
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Consider a dilute monoenergetic particle beam that propagates in a dense cold plasma. The rest mass of each beam 
particle is mb and its energy is mb7, as measured in the background plasma frame0 In this frame the beam momenta 
are narrowly collimated within a half opening angle ~ F"^, where 7 ^ F ^ 1. The following calculation makes no 
assumptions about the beam and plasma composition. 

Since we consider a monoenergetic beam, we can neglect the longitudinal momentum spread of the beam, as seen 
in the upstream frame, in which the calculations will be conducted. The transverse momentum spread in this frame 
is ~ TOb7/r. In order to investigate the stability of such a beam, we consider a system containing a beam of density 
rib and mean velocity Vh ^ c that interpenetrates a cold background plasma of density ripi ^ Tib. Both densities 
are measured in the upstream plasma rest frame (in the shock frame the beam density is smaller by a factor of F 
while the upstream plasma density is larger by the same factor). The beam density accounts for all species (e.g., in 
a proton-electron beam, we have Uh = n\^^e + ?^b,p), while the plasma density includes only the lightest species (e.g., 
in a proton-electron plasma, we have Upi — Tipi ^ and in a pair plasma, we have ripi — Jipj e- -I- '^pi,e+)) since the 
heavier species is assumed to be static. The system is charge neutral so plasma drift neutralizes the beam charge. In 
a proton-electron plasma and beam, the electrons drift with velocity Vpi ~ C(nb/?ipi)fb and cancel the beam current, 
where 

~ nh 

is the fractional beam charge satisfying — 1 < C < 1 . We show later that this drift has a negligible effect on the 
instability spectrum and on the growth rate. 

Given the symmetry of the dispersion relatio n with respect to t he rotation around the parallel direction, a two- 
dimensional survey of the fc-space is required (|Bret. Firpo fc Deut sch 2005). In general, such calculation requires 
solving the linearized Maxwell- Vlasov equations. The electrostatic, or longitudinal, approximation, i.e., k x E = 0, is 
not appropriate in general because the modes evolve continuously from k x E = along the beam direction, to k x E 7^ 
in the perpendicular direc tion. For example, the filamenta tion instability is found for k _L Vb and has k x E 0, but 
it need not have k • E = ()Bret. Gremillet fc Bellidoll2007D . Thus, in general, the electrostatic approximation can lead 
to a neglect of the most unstable modes. However, in the specific case of interest, where t he be am is relativistic and 
dilute, oblique or parallel parts of th e spectrum govern the linear phase of the instabilit y (jBret al. I 2008) and the 
electrostatic approximation is valid (jBret. Firpo fc DeutschI [200I : iGremillet et al. I [20061) . Thus, in what follows, we 
adopt the electrostatic approximation. 

The longitudinal dispersion equation thus reads ([Ichimarul 1 1 9 73[) . 

e(k,.)^l + E^ f^^^^d^P = 0, (2) 

where the sum runs over the species a which are here the beam species, plus the lightest species of the cold background 
plasma. We now choose the z axis as the direction of beam propagation, and label x the normal direction. The 
momentum distribution function of each beam species reads: 

/b(p) = ^^{P^ - 7"ib«b) + Pt) - Q{p. - Pt)] , (3) 

where Pt = ^m^Vh/T, while the distribution function of each (background) plasma species, is 

/pi(p) = KPz + Cmp\Vp\)5{p^), (4) 

where m-p\ is the mass of the lightest plasma particle, typically an electron, and Q{t) is the Heaviside step function. 
Note that 7TOb is equal for all beam species and that in a pair plasma C — 0. 
The dispersion equation obtained when inserting these functions in equation ^ can be written concisely as 

(u — CZzay mb 

where 

^ ^ _Jl ^ z,^ ~ ptj\z, - ptz,) 

Zl + Zl[Z,-^cj){p)f -Zlp-^ 2i^{Zl + Zl) J_^ct^{pt)[Z, + ptZ,-Cjct^{pt)f 
Here 7 is the Lorentz factor of beam particles with mass mb and we define the following symbols 

?T-b „ kc 1 W , , . 1 , 77-, ^ , , 

a = — , Z = — , p=-, Cj=—, <^(t) = -v/1 + 7' 1 + t2 , 7 
n-pi Wpi r Wpi 7 

where Wpj = Airupie'^ / rric is the plasma frequency of the background lightest particles and we approximate Vh w c. The 
second term on the r.h.s. of equation ([SJ pertains to the lightest particle species of the background plasma, and the 
last term accounts for the beam response. 

^ If the beam contains several particle species, then we assume energy equipartition between particles of different masses. For example, 
in a proton-electron beam rrae7b,e = "*p7b,p (= ™b7)- Given the highly relativistic setting considered here, the slight velocity drift between 
species of different mass (e.g., since 7b, e 7^ 7b, p) can be neglected. 
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Fig. 1. — Growth rate in units of the plasma frequency as a funetion of the dimensionless wave vector Z = kc/ojpi, in the case of a pure 
electron beam. Parameters are a = 10-2, 7 = 100, and T = 5. 




0=tan~''(k^/k^P 

Fig. 2. — Growth rate in units of the plasma frequency as a function of the mode angle tan"^ (fcx/fc|| ). Parameters are the same as in 
Figure [T] 



Calculations can be simplified in the present regime noting that < p 1 and 7^1, which implies that (/'(p) ~ 1 
and 4'{pt) ^ 1 for — 1 < t < 1. The quadrature in equation ([5]) can be further simplified to obtain 



All types of beams and plasmas exhibit a similar unstable spectrum. The reason is that the only dimensionless 
parameter that cannot be removed through scaling is the fractional beam charge C . It is the beam charge density that 
determines whether the background plasma drifts and what is the direction of the drift. However, the drift velocity is 
proportional to the density ratio a and is thus very small. Therefore in equation ([S]), the parameter C is only found 
in a product with the density ratio a and the effect of space charge in the beam is negligible (note that the value of 
the third factor in the product \s ^ \ for the most unstable mode). 

A typical two-dimensional map of the unstable spectrum for a pure electron beam is shown in Figure[Tl for a ~ 10^^, 
7 = 100, and F = 5. The growth rate of the fastest growing unstable mode as a function of the mode angle 
tan^^(fcx/fc||) is shown in Figure [5] The spectrum exhibits two local maxima. One corresponds to the parallel mode 
with {Z\\^Z\) sd (1,0) and the other to the oblique mode with (Z^^Z\_) « (1.2,1). The relative height of the two 
maxima, and hence the character of the fastest growing mode — whether it is parallel or oblique — depends on the 
specific values of the parameters. Also noticeable is the narrowness of the range of wave vectors containing unstable 
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modes. With a high beam Lorentz factor, a mode can only grow out of its interaction with the beam if a certain 
resonance condition is strictly fulfilled (we discuss this condition in ii l3.2p . Modes which are "out of tune" are quickly 
damped away. 

The total absence of unstable modes in the perpendicular direction is not always just a direct consequence of 
the electrostatic approximation, which does not account for such non-longitudinal modes. The theory of unstable 
filamentation modes for the present setting shows that they are completely shut down if the transverse beam angular 
spr ead as seen in the upstream frame ^ is larger than \amp] / (^mY,)V^^ , that is, when the density ratio is limited 
by ()Silva et al. Il2002t iBret. Firpo fc DeutschlbOOa iRabinak. Katz fc WaxmanllMO) 

a < — (filamentation suppression). (9) 

Since we work here in the regime in which a ^ 1 <C F ^ 7, filamentation modes are unstable only in a small portion of 
the phase space, especially in an ion-electron plasma. Moreover, if the first shock-crossing cycle of the DSA accelerates 
particles to Lorentz factors 7 ^ (e.g.. Gallant & Achtcrbcrg 1999), then the filamentation modes are stable in the 
entire precursor for any type of plasma. Finally, let us emphasize that even when filamentation modes are unstable, 
their growth rate is slower than the fastest growing two-stream mode that we find here, so they do not dominate 
anyway. 

The growth rate of the fastest purely parallel growing mode can be found analytically. The dispersion equation for 
such wave vectors, Z — Z^z, is easily derived from equations ([2]) and ^ and reads 

1 TOoi (3 + 7^^) , , 

= 1 a—— ' (10) 

{u- CZ.af mb 3(tD - Z.^-^^ ^ ' 

This relation is mathematically similar to the dispersion relation of the two-stream instability. The maximum parallel 
growth rate Jpar is reached for ~ 1 and reads 

3I/6 J a 

The growth rate of the fastest growing oblique mode can be solved for numerically. We conducted a systematic 
numerical search for the fastest growing mode for representative choices of a, 7, and F. We find that when a 1 and 
1 ^ F ^ 7, the maximum oblique growth rate Sohi reads 

^obi ~ 7; — '^pi- (12) 

3 Vm-b 7 / 

Equating the two growth rates yields the criterion for the dominance of each, so that the global fastest growing 
mode grows at the rate 
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F > 1. 



(13) 



Some of the general trends evident in equations (fTTj) and (fT2|) were to be expected. Namely, the growth rate weakens 
when the beam density decreases or when the inertia of the beam particles increases. The fact that the growth rate has 
two local maxima with opposite dependencies on F is nontrivial and is discussed in th e next section. It is interesting to 
compare this result to the scenario in which a cold beam interacts with a cold plasma (jFamberg. Shapiro fc Shevchenkol 
119701 ). There, the fastest growing mode is always oblique and its growth rate scales as (a/7)^/'^. 

Finally, we stress that equations (jlip and (|12p provide the growth rates for monoenergetic beam with an idealized 
distribution function given in equation (jS]). The true distribution function in collisionless shocks with DSA is of course 
more complicated. Since the exact value of the growth rate and its dependence on the various parameters depend 
on the exact shape of the distribution function, we expect the result to be somewhat different if a more realistic 
distribution is considered. However, we still expect that the key proportionalities evident in equations ()11|) and (|12[) 
persist in the general case, although the exact values of the power-law indices may vary. We also expect that there 
will still be a competition between the parallel and oblique modes. 



3.2. Intuitive physical description 

The oblique mode is a two-stream-like electrostatic mode, and as such it can be understood with the aid of the 
intuition gained from the consideration of a one-dimensional, longitudinal two-stream mode. In one dimension, a wave 
mode interacts with near-resonant particles that move with velocities similar to the phase velocity of the mode. These 
particles "surf the wave in the sense that they spend a long time in phase with the wave. During surfing, the velocity 
of near-resonant particles further converges to the phase velocity of the wave. The particles that are slightly faster 
than the wave tend to transfer energy to the wave, and thus excite it, while those that are slightly slower than the wave 
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Fig. 3. — Illustration of the beam and wave mode geometry. 
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Fig. 4. — Growth rate map in the plane defined by ^"^"^ [cos (6* + F ^)] ^ and 6 for unstable modes with < Z± < 2. The color 
shows the growth rate. Parameters are a = 10-2, 7 = 1000, and F = 100. The narrow concentration of the unstable modes around the 
resonant condition ^"^"^ [cos(g + F" -'^)]--' = 1 is clearly evident. 



tend to receive energy from it, and thus damp the wave. Since growing modes are those with a net energy gain, a mode 
is unstable if the number of particles moving slightly faster than the wave exceeds that of those moving slightly slower. 
Thus a distribution function f{v) should be two-stream unstable only if (df /d\v\)^ ^ > 0, where Uphaso — Re(a;)/fc 
is the phase velocity of the wave mode. 

The three dimensional case is similar, except that now it is Wk, the component of the velocity parallel to the wave 
vector k, that determines the stability of the mode. In the configuration under consideration all the accelerated 
particles move at about the speed of light and have momenta that are coUimated to within the small angle ^ from 
the direction of shock propagation (see Figure |3|). For any generic wave vector, except for the almost perpendicular 
ones, the range of particle velocities projected onto the wave vector does not include the zero velocity. Thus a mode 
with phase velocity equal to the minimum of the projected beam particle velocities has many particles that are moving 
slightly faster than the wave but none that are moving slower than the wave. Based on the one-dimensional analysis 
in the preceding paragraph, the mode is unstable. Neighboring modes with an excess of particles moving faster than 
the wave should also be unstable. All unstable two-stream oblique modes satisfy the resonance criterion 



Wpha 



Re(w) 

ck 



(14) 



where 9 is the angle between k and the z axis. This include the parallel mode, 6 = 0, and oblique modes up to 

6 = tt/2 — r^^. Figure m shows a growth rate map in the plane of ^^^^^[cos(0 -I- F"^)]^^ and 9, where the resonance 
is clearly evident. The resonance condition for instability in the regime that we consider, up to terms of the order of 
r~^, can therefore be written as 

Re(w) « ckn - (15) 



8 



TWO-STREAM INSTABILITY IN RELATIVISTIC SHOCKS 



This intuitive description does not select the fastest growing mode among all the unstable modes. In ti3.1l we analyzed 
the dispersion relation to pin down the dominant mode and derived the dependence of growth rate on the shock Lorentz 
factor r in the parallel and oblique regimes. The different dependencies of the growth rates on F can be understood 
as follows. In general, the growth rate depends on two competing and opposite factors. One is the number of particles 
that are in resonance with the mode. The energy of these particles can in principle be tapped toward contributing 
to the growth of the mode. The other relates to efficiency with which energy transfer from the particles to the mode 
can take place. The second factor is the reason that oblique modes play an important role in relativistic beams. It is 
much easier the reduce the projected velocity of an ultra-relativistic particle in the direction of k, and thus to trap the 
particle in a resonance with the mode and transfer its energy to the mode, by deflecting the particle sideways, than 
by decelerating the particle. From this it follows that larger value of the angle between the resonant particles and the 
wave, 9 + T~^ , contribute to increasing the growth rates. On the other hand, now returning to the first aforementioned 
factor, larger values of F and smaller values of 9 imply a distribution function that is more compact when projected 
onto the direction of k. This ensures that a larger number of nearly- resonant particles is available. From this, however, 
it follows that having a larger value of 9 + T^^ also works in the opposite direction, namely, toward decreasing the 
growth rates. It is the balance between these to opposite dependencies of the two effects on 6* -I- F~^ that determines 
the location in k space of the fastest growing mode. Now, in oblique modes, where 9 3> F~^ and thus 9 + F~^ Ki 9, 
increasing F only increases the compactness of the distribution function in the direction of k, thereby increasing the 
growth rate of the fastest growing oblique mode. In parallel modes where w and thus 9 + F~^ w F~^, increasing 
F increases the number of resonant beam particles but makes it harder to exploit their energy. It turns out the latter 
factor wins and the growth rate increases when F is smaller. 

4. APPLICATION TO EXTERNAL SHOCKS IN GRBS 

Here we assess under which conditions the two-stream-like instability in the precursor of GRB external shocks 
can become nonlinear and strongly affect the shock structure, assuming that the acceleration process is DSA. We 
do so by comparing igrow in such a shock (see the definition in §[2|) to l/5max found in § 13.11 A GRB external 
shock is a spherical relativistic collisionless shock that propagates in a weakly magnetized pro ton-electron plasma . 
Electron ac celeration in GRB external shocks is expected to be limited by the radiative cooling (|Li fc WaxmanI [20061 : 
iCouch. Mil osavlievic & Nakar 2008) and thus the precursor beam should contain mostly protons at higher energies. 
The accelerated particle spectrum is taken to be a power-law dNi^cc/dj oc where A^acc is the total number of 
accelerated particles in the shocked region and the power law index is in the range 2 < p < 3. We assume that the 
minimum accelerated particle energy is ^ F^, where F is the shock Lorentz factor, as expected for the accelerated 
particles that have completed their first full Fermi cycle. The density of the accelerated particle precursor can then be 
parameterized by the fraction Cacc of the shock energy that goes into the accelerated particles. The total accelerated 
particle number is then A'acc ~ ^accE / (T^mpC^), where E is the total shock energy. The density of accelerated particles 
with Lorentz factors ~ 7, as measured in the upstream frame, is nacc(7) ~ (r^/i?^)(7/F^)^~P7Vacc, where R is the shock 
radius. We assu me that the density n of th e unshocked upstream is constant. Hydrodynamics of an ultra-relativistic 
blast wave fe.g.. lBlandford fc McKejll976( ) relates n to the blast wave energy E, radius R, and Lorentz factor F; with 
this, one finds that the density ratio measured in the upstream frame, a = ng^cc/n^ is 

Note that a retains no explicit dependence on E or n. At any distance from the shock the accelerated particle density 
is dominated by the lowest energy particles that reach that distance. 

We first determine which among the oblique and parallel modes dominates. Substituting equation (1161) into (1131) . 
using rrib = "mp, we find that with 7 ~ 10^ (the dependence on the particle Lorentz factor is weak), for F > 60 the 
oblique mode dominates while for F <C 60 the parallel mode dominates. Since most of the afterglow observations cover 
the evolution when F < 60, in what follows, we will focus on the low-F regime in which the parallel mode dominates. 
The growth rate derived in equation (ITT|) decreases with increasing 7 and scales as (5max oc (a/7F^)^/^ oc j~p/^r'^(p~^y^ . 
The time available for growth is tgrow ^ ^{l)/c where A(7) is the distance from the shock that particles with Lorentz 
factors 7 can reach. The maximum distance is Amax ~ i?/(10F^), and it is reached by the protons that are accelerated 
to the highest energies 7max- A pre-existing shock upstream magnetic field with a strength of a few ^G and with a 
coherence length that is larger than the shock radius, as expected in the ISM, can accelerate protons up to 7max ^ 10^. 
If accelerated particle streaming in the shock precursor can somehow amplify the pre-existing field to the point of 
equiparti tion with the rest energy of the circum-burst medium, proton acceleration to energies 7max ~ 10^° may be 
possible (Milos avlievic fc Nakai][2006l ) . We approximate the dependence of the distance from the shock on the particle 
Lorentz factor with a power-law A(7) oc 7". If all the accelerated particles are defiected back into the shock by a 
magnetic field with coherence length larger than R, then we have s = 1. If, on the other hand, they are defiected by 
a small-scale with coherence length As ^ R, then we have s = 2 at all distances A > As. 

The maximum number of e-foldings by which an unstable mode can grow is the ratio igrow/'^maxj implying that 
A^c-foid(7) = '5niax(7)A(7)/c. Using equations (HI]) and ([T6|) . we find 
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The dependence on E and n is rather weak, while the scahng with varies from being hnear to being quadratic 
for 2 < p < 3. This imphes that if the accelerated protons carry a non-negligible fraction of the total energy of the 
shock, e.g., eacc ^ 10"'^, but the maximum proton energy is relatively low, 7„iax ^ 10^, then as long as the external 
shock remains relativistic A''e-foid(7max) 3> 1, and at the largest distance from the shock there is enough time for the 
two-stream modes to become nonlinear, then the instability can strongly influence the state of the upstream plasma 
and accelerated particle beam dynamics. If, in the other extreme, 7rnax ~ 10"'^'' and s = 1, then -/Vo-foid(7max) ^ 1 
throughout the external shock, and the two-stream instability should not affect its structure. For intermediate values 
of 7max or when s = 2, the number of e-foldings is below unity when F is large, but it can become larger than unity 
as F decreases; this transition may imply a change in the shock structure. Finally, we note that for low values of 7, 
i.e., close to the shock, and for high values F, it is possible that the accelerated particle density exceeds the upstream 
plasma density, a > 1, in which case our analysis is no longer valid. 

5. SUMMARY 

We explored the linear growth rate of electrostatic two-stream-like modes that grow when a dilute ultra-rclativistic 
beam, which is relativistically hot in its own rest frame, propagates in a cold, unmagnetized or weakly magnetized, 
background plasma. This beam-plasma configuration describes the precursor containing accelerated particles that runs 
ahead of an ultra-relativistic, unmagnetized or weakly magnetized, coUisionless shock. We considered various types of 
plasma (electrons and ions or pairs) and beams of various compositions and studied the fastest growing kinetic modes. 
Such modes may thus be the principal agent of nonlinear pre cursor- upstream int eraction that is expected in these 
shocks and that is also apparent in numerical simulations (e.g.. iKeshet et a7~ll2009f) . 

Our main finding is that for any beam and plasma composition and for any shock and beam parameters, as long as 
7 3> F ^ 1 and the background plasma is cold, there exist two-stream-like unstable modes. These are the modes that 
satisfy the resonance criterion Re(a;) /ck cos{9 + F"^), where 6 is the angle between wave vector k and the direction 
of beam propagation. This includes the parallel mode with 9 = 0, and oblique modes up to 6* = 7r/2 — F~^. The 
resonance criterion can also be written as Re(a;) w ck^^ — ck±/T. This behavior departs from that of the filamentation 
(or transverse Weibel) modes, which are suppressed over a large range of the parameter space. Moreover, in the 
interaction of a dilute accelerated particle precursor with a dense, cold, unmagnetized, and unperturbed upstream, the 
fastest growing electrostatic two-stream-like mode always grows faster than any of the filamentation modes when the 
latter are unstable. We therefore do not expect filamentation modes to play a role in the shock precursor containing 
accelerated particles, but possibly only in the shock transition itself. 

We find that the spectrum of unstable electrostatic modes exhibits two local maxima. One contains parallel modes, 
with wave vectors (fcy ,k±)^ —(1, 0), and the other contains oblique ones with (fcy ,k±)^ ^(1+F^^, 1). The relative 
height of the two maxima, antf therefore the character of the dominant mode — whether it is of the parallel or oblique 
type — depends on the parameters and is most sensitive to F. For large values of F, an oblique mode dominates, while 
for small values of F, it is the parallel mode that dominates. The criterion for finding the dominant mode exhibits 
a similar trend with TOpia/(mb7), but this dependence is much weaker. Nevertheless, in pair or electron beams, the 
oblique mode dominates over a larger portion of the parameter space than it does in ion beams. 

The time that is available to an unstable mode for its growth is limited in the precursor of astrophysical ultra- 
relativistic shocks. Nevertheless, if sufficient time is available for the fastest growing mode to become nonlinear, the 
nonlinear outcome of the instability may significantly affect the shock structure. It can lead to momentum transfer 
from the accelerated particles to the shock upstream, thereby modifying the hydrodynamic structure of the shock. We 
compared the growth time of the fastest growing two-stream-like mode to the time available for its growth in GRB 
external shocks. We find that for all except for the earliest phases, when the shock Lorentz factor has dropped below 
~ 100, the parallel mode dominates. If the maximum Lorentz factor of the accelerated protons is modest, 7max ~ 10^, 
then this mode has enough time to become nonlinear, while if 7max is significantly larger, then the mode may become 
nonlinear only for relatively low values of F. We conclude that the electrostatic two-stream-like instability may play 
an important role in accelerating particles and in shaping of the structure of the external shocks in GRBs. 
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